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Abstract This review summarizes the ohservational constraints on the mass spectrum of 
compact objects in X-ray binaries. We currently have 20 X-ray binaries with 
confirmed black holes, based on dynamical information (i.e. mass in excess of 3 
Mq). In two cases, V404 Cyg and GRS 1915-1-105, the black hole mass exceeds 
the maximum predicted by current Type Ib supemovae models and challenges 
black hole formation scenarios. The great majority of black hole binaries are 
members of the class of X-ray Transients, where long periods of quiescence en¬ 
able spectroscopic studies of the faint donor stars. On the other hand, neutron 
star binaries are mostly found in persistent binaries, where reprocessed light 
from the accretion disc overwhelms the companion star and precludes mass es¬ 
timates. New results, based on the detection of optical fluorescent lines from 
the donor star and X-ray burst oscillations, provide the best prospects for mass 
constraints of neutron stars in persistent X-ray binaries. 
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Introduction 

Building the mass distribution of compact objects is a fundamental experi¬ 
ment in modem Astrophysics which can only be done in X-ray binaries. The 
interest of this research is two-fold: set constraints on the equation of state of 
nuclear matter (hereafter EOS) and test models of supernovae explosions and 
close binary evolution. X-ray binaries are interacting binaries where a normal 
star transfers matter onto a compact object, a black hole (BH) or neutron star 
(NS). Matter is accelerated in the strong gravitational field of the compact star 
and heated up to ~ 10^ K before being accreted. This is the canonical model, 
first proposed by Shklovskii (1967) to explain the new X-ray sources detected 
in the 60’s and 70’s by X-ray satellites such as UHURU or Einstein. In the 
80’s and 90’s, a new generation of higher sensitivity satellites, such as Ginga, 
Rosat and XTE, revealed a population of several hundred X-ray binaries in 
the Galaxy with X-ray luminosities in the range 10^® - 10^® ergs s“^ and dis¬ 
playing rapid variability, down to the kilohertz regime. And nowadays, the 
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large collecting area and high angular resolution of satellites such as Chandra 
and XMM-Newton allow us to resolve the X-ray binary population in nearby 
galaxies, obtaining the X-ray luminosity function and proving that they are 
responsible for typically a third of the galaxy’s total X-ray emission. 

X-ray binaries are classically divided into two populations, based on the 
nature of the optical (companion) stars: high mass X-ray binaries (HMXBs) 
and low mass X-ray binaries (LMXBs). In HMXBs the companion is a hot 
luminous 0-B supergiant whose bolometric luminosity completely dominates 
the energy spectrum. These binaries are relatively young (short lived) with 
estimated ages in the range 10^ — 10® yrs. There are about ~ 100 such binaries 
in the Galaxy, strongly concentrated towards the spiral arms and, hence, they 
are considered as good tracers of the star formation rate (Grimm, Gilfanov & 
Sunyaev 2002). Neutron stars in HMXBs are recognized by the presence of 
regular X-ray pulsations as a consequence of their strong magnetic fields. 

On the other hand, LMXBs contain low-mass companion stars of spectral 
types (mainly) K-M. These binaries are very compact, with orbital periods 
clustering around 4-8 hours and concentrated in location towards the galac¬ 
tic bulge. They are associated with the old stellar population and show a 
large spread in galactic latitude, interpreted as a signature of kick velocities 
received when the supernova explosion formed the compact star. Neutron stars 
in LMXBs are also revealed by the exhibition of sporadic X-ray bursts i.e. 
thermonuclear eruptions of matter slowly accreted over the surface. There 
are ~ 200 X-ray active or "persistent" sources in the Galaxy, and about 10® 
"transients", which only show X-ray activity occasionally. The reason for this 
transient behaviour is the mass transfer rate from the donor star M 2 , which is 
driven by the binary/donor evolution. In X-ray transients, M 2 is lower than a 
critical value ~ 10 “® M 0 yr“® and this triggers thermal-viscous instabilities in 
the accretion disc. This causes enhanced mass transfer episodes onto the com¬ 
pact object, the so-called "outbursts", with recurrence times of a few decades ( 
King 1999). In the interim, these systems remain in the "quiescent" state, with 
typical X-ray luminosities below ~ 10®^ ergs s“®. It is during these periods 
when we can attempt to detect the faint low-mass star and extract dynami¬ 
cal information. The many studies performed on quiescent X-ray transients 
have demonstrated that BHs outnumber NS by more than 70%, and hence the 
transients provide excellent hunting grounds for BHs. More details about the 
optical and X-ray properties of galactic X-ray binaries can be found in Charles 
& Coe (2004), McClintock & Remillard (2004) and Psaltis (2004). 

1. Establishing Black Holes 

As opposed to NS, which often show X-ray bursts or pulses, the best ob¬ 
servational evidence for a BH is still its mass function /(Mx). This equation 
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relates the masses of the eompaet objeet M^, the eompanion star Me (or alter¬ 
natively the binary mass ratio q = Mc/Mx) and the binary inelination angle i 
with two quantities to be extraeted direetly from the radial veloeity eurve of the 
eompanion star: the orbital period Porb and the radial veloeity semi-amplitude 
K: 


\ _ Ml sin^ i _ Mx sin^ i 

27rG “ (Mx + M,f “ (l + g)2 

Sinee Me > 0 and 0 < i < 90° it is straightforward to show that f{M) 
is a lower limit on M^. Therefore, a mass funetion larger than 3 M©, i.e. the 
maximum gravitational mass of a NS (Rhoades & Ruffini 1974), is taken as a 
seeure proof for the existenee of a BH, independently of the aetual values of i 
and q. 

Table 1 presents an updated list of eonfirmed BHs based on this simple dy- 
n ami cal argument. We currently have 20 BHs, with orbital periods ranging 
from 4.1 hours to 33.5 days. There are 17 transient LMXBs and only 3 per¬ 
sistent systems, the HMXBs Cyg X-1 and two sources from the Large Mag¬ 
ellanic Clouds: LMC X-1 and LMC X-3. Note that, in addition to the three 
HMXBs, six transient LMXBs have mass functions < 3 Mq. However, we 
have solid constraints on the inclination and the companion’s mass for these 
binaries which result in Mx > 3 M©. The last one added to the list, BW Cir, 
contains an evolved G5 donor star in a 2.5 day orbit. Its optical luminosity 
places the binary at a distance > 27 kpc and makes BW Cir the furthest BH 
binary in the Galaxy yet. Remarkably, its large systemic velocity (103 km s“^) 
is in good agreement with the projected velocity of the Galactic differential 
rotation at that distance (see more details in Casares et al. 2004a). 

In addition to the mass function, one obviously needs a knowledge of the 
inclination and the mass ratio to fully determine the masses of the two stars. In 
the case of a BH binary, we face a "single-line spectroscopic binary" problem 
and all the information must be extracted from the optical companion. How¬ 
ever, there is a complete solution to the problem which involves: (i) the deter¬ 
mination of the binary mass ratio through the rotational broadening Vrot sin i 
of the companion’s absorption lines and (ii) the determination of the inclination 
angle by fitting synthetic models to ellipsoidal lightcurves. This is the classic 
method to derive masses. Further details of these techniques and the systemat- 
ics involved can be found in several review papers e.g. Casares (2001), Charles 
& Coe (2004). 

Following this prescription, we currently have 15 reliable BH masses which 
are listed in the last column of Table 1, updated from the compilations in Grosz 
(2003) and Charles & Coe (2004). Figure 1 plots these masses with their la 
errorbars. BH masses spread between 4 and 14 Mq, with typical uncertainties 
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Table 1. Confirmed black holes and mass determinations 


System 

f* orb 

(days) 

f(M) 

(Mq) 

Donor 
Sped. Type 

Classification 

/Hx 

(Mq) 

GRS 1915+105 

33.5 

9.5 ± 3.0 

K/MIII 

LMXB/Transient 

14 ±4 

V404 Cyg 

6.470 

6.08 ± 0.06 

KOIV 

,, 

12 ±2 

Cyg X-1 

5.600 

0.244 ± 0.005 

09.7 lab 

HMXB/Persistent 

10 ±3 

LMC X-1 

4.229 

0.14 ±0.05 

07 III 

,, 

>4 

XTEJl 819-254 

2.816 

3.13 ±0.13 

B9III 

LMXB/Transient 

7.1 ± 0.3 

GRO J1655-40 

2.620 

2.73 ± 0.09 

F3/5 IV 

,, 

6.3 ± 0.3 

BW Cir “ 

2.545 

5.75 ± 0.30 

G5IV 

,, 

> 7.8 

GX 339-4 

1.754 

5.8 ±0.5 

- 

,, 


LMC X-3 

1.704 

2.3 ± 0.3 

B3 V 

HMXB/Persistent 

7.6 ± 1.3 

XTE J1550-564 

1.542 

6.86 ±0.71 

G8/K8 IV 

LMXB/Transient 

9.6 ± 1.2 

4U 1543-475 

1.125 

0.25 ± 0.01 

A2 V 

„ 

9.4 ± I.O 

H1705-250 

0.520 

4.86 ±0.13 

K3/7 V 

,, 

6±2 

GS 1124-684 

0.433 

3.01 ±0.15 

K3/5 V 

,, 

7.0 ± 0.6 

XTE J1859+226 ** 

0.382 

7.4 ± 1.1 

- 

,, 


GS2000+250 

0.345 

5.01 ±0.12 

K3/7 V 

„ 

7.5 ± 0.3 

A0620-003 

0.325 

2.72 ± 0.06 

K4 V 

,, 

II ± 2 

XTE J1650-500 

0.321 

2.73 ± 0.56 

K4 V 

,, 


GRS 1009-45 

0.283 

3.17 ±0.12 

K7/M0 V 

„ 

5.2 ± 0.6 

GRO J0422+32 

0.212 

1.19 ±0.02 

M2 V 

,, 

4± I 

XTE Jill8+480 

0.171 

6.3 ± 0.2 

K5/M0 V 


6.8 ± 0.4 


“ The 1-year alias period at 2.564 days is equally significant. In this case the BH would be even strengthen 
with /(Mx) = 6.60 zb 0.36 Mq (see Casares et al. 2004a). 

^ Period is uncertain, with another possibility at 0.319 days (see Zurita et al. 2002). This would drop the 
mass function to /(Mx) = 6.18 Mq. 

^ Mass function updated after Torres et al. (2004). 


in the range 5-33%. Also in Fig. 1 we show 33 well determined NS masses, 
extracted from Stairs (2004) and Lattimer & Prakash (2004), also with la un¬ 
certainties. The most precise NS masses have been measured in a group of bi¬ 
nary radio-pulsars. They are descendants of HMXBs, composed of two young 
pulsars whose orbits are known to great accuracy from pulse time delays. Rela¬ 
tivistic effects lead to NS mass determinations with exquisite accuracy and they 
display a normal distribution centered at the canonical value of 1.35 Mq with a 
very small dispersion of ±0.04 Mq. Dynamical masses are also available from 
pulsing NS in seven HMXBs, six of which are eclipsing. However, the uncer¬ 
tainties are much larger because the radial velocities of the optical stars are 
distorted by non-Keplerian perturbations, caused by their strong winds. A few 
measurements also exist for LMXBs and binary millisecond pulsars (BMPs). 
BMPs are descendants of LMXBs, composed of a millisecond NS (spun up 
by accretion) and a detached white dwarf. Both LMXBs and BMPs are poten- 
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Figure 1. Mass distribution of compact objects in X-ray binaries. The arrows indicate lower 
limits to BH masses. The dotted vertical line marks the maximum mass allowed for a stable NS. 


tial sites for massive NS beeause binary evolution prediets aeeretion of several 
tenths of a solar mass during their lifetimes. 

A fundamental result for the understanding of nuelear matter would be to 
find a NS more massive than 1.6 Mq since this would rule out soft EOS ( 
Brown & Bethe 1994). Currently the best candidates are found in the LMXB 
Cyg X-2 (1.78 ± 0.23 Mq: Casares, Charles & Kuulkers 1998, Orosz & Ku- 
ulkers 1999), the BMP J0751-bl807 (2.2 ± 0.2 M©: Nice, Splaver & Stairs 
2004) and the HMXBs 4U 1700-37 (2.44 ± 0.27 M©: Clark et al. 2002) and 
Vela X-1 (1.86 ± 0.16 M©: Barziv et al. 2001). The latter case, however, 
is less secure since the radial velocity curve of the companion is affected by 
systematic excursions which prevents a confirmation of the mass estimate (see 
Barziv et al. 2001). In the case of J0751-I-1807 and 4U 1700-37, soft EOS are 
ruled out even at the 95% level. 

Very recently, large NS masses have also been reported for 2S0921-630 
(=V395 Car) by two different groups: 2.0-4.3 M© (Shahbaz et al. 2004) and 
1.9-2.9 M© (Jonker et al. 2004). However, the compact object in this EMXB 
has never shown any evidence for X-ray bursts or pulsations, so it could be a 
low-mass BH. This is also the situation for the compact object in 4U 1700-37, 
since it has never shown any NS signature. 
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Compact Object Mass (M^) 


Figure 2. Observed mass distribution of compact objects in X-ray binaries (shaded his¬ 
togram), compared to the theoretical distribution computed in Fryer & Kalogera (2001) for 
"Case C -I- Winds" scenario (dotted line). The model distribution has been re-scaled for clarity. 


Figure 2 presents the histogram of eompaet objeet masses eompared to the 
theoretieal distribution of remnants eomputed in Fryer & Kalogera (2001) for 
the ease of binary interaetion under Case C mass transfer (i.e. Common Enve¬ 
lope evolution after eore helium ignition) and mass loss through winds in the 
Wolf-Rayet phase. This is the most realistie seenario sinee evolution through 
Case B mass transfer (i.e. Common Envelope and H envelope removal before 
eore helium ignition) fails to produee BH remnants > 3 M© with the eurrent 
Wolf-Rayet mass loss rates (Woosley, Eanger & Weaver 1995). The model 
prediets a BH mass eut at 12 Mq whieh is diflieult to reeoneile with the high 
masses measured in V404 Cyg (12 ± 2 Mq) and GRS 1915-1-105 (14 ± 4 Mq). 
Despite the large uneertainties in the masses, these two X-ray binaries seem to 
pose a ehallenge to BH formation theories and, in partieular, suggest that mass 
loss rates in the WR phase are overestimated. 

Our histogram also shows a shortage of objeets at 3-4 Mq, not predieted by 
the model distribution. If the gap is eventually eonlirmed, it eould strongly re- 
striet the supernova explosion energy sinee, as explained in Eryer & Kalogera 
(2001), it ean be reprodueed by a step funetion dependenee with the progeni¬ 
tor’s mass. However, seleetion effeets eould also be playing a role here sinee 
low mass BHs are likely to show up as persistent X-ray sourees, where dynam- 
ieal masses are diffieult to obtain. The ease of 2S0921-630 eould well be an 
example. Obviously we are limited by low number statisties in the observed 
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distribution of compact remnants. Clearly more X-ray transient discoveries 
and lower uncertainties in the mass determinations are required before these 
issues can be addressed and the form of the distribution can be used to con¬ 
strain supernova models and X-ray binary evolution. 


2. Mass Determination in Persistent LMXBs 

So far we have been dealing with mass determination in quiescent X-ray 
binaries. Dynamical studies in persistent LMXBs are, on the other hand, ham¬ 
pered by the huge optical luminosity of the accretion disc. This is driven by 
reprocessing of the powerful (Eddington limited) X-ray luminosity and com¬ 
pletely swamps the spectroscopic features of the faint companion stars. 

New prospects for mass determination have been opened by the discovery of 
high-excitation emission lines arising from the donor star in Sco X-1 ( Steeghs 
& Casares 2002). The most prominent are found in the core of the Bowen 
blend, namely the triplets NIII A4634-40 and CIII A4647-50. In particular, the 
NIII lines are powered by fluorescence resonance which requires seed photons 
of Hell Ly-a. These narrow components move in phase with each other and 
are not resolved (FWHM=50 km s“^ i.e. the instrumental resolution), an in¬ 
dication that the reprocessing region is very localized (Fig. 3). The extreme 
narrowness rules out the accretion flow or the hot spot and points to the com¬ 
panion star as the reprocessing site. 



Figure 3. Trailed spectra of the nar¬ 
row CIII+NIII Bowen emission lines and 
Hell A4686 in Sco X-1. After Steeghs & 
Casares (2002). 


Figure 4. Radial velocity curve of the 
sharp CIII+NIII Bowen lines (top) and the 
wings of Hell A4686 (bottom) in Sco X-1. 


The radial velocity curve of the donor star can be extracted through a com¬ 
bined multigaussian fit to the three CIII/NIII lines (top panel of Fig. 4) or 
using more sophisticated Doppler tomography techniques (see e.g. Casares et 






al. 2003). Furthermore, the veloeities are in antiphase with the wings of the 
Hell A4686 emission, whieh approximately traee the motion of the eompaet 
star (bottom panel of Fig. 4). This work represents the first deteetion of the 
eompanion star in Seo X-1 and opens a new window for extraeting dynami¬ 
cal information and deriving mass functions in the population of ~20 LMXBs 
with established optical counterparts. 

Follow-up campaigns, using the AAT, NTT and VLT telescopes, have en¬ 
abled us to extend this analysis to other fainter LMXBs leading to the detec¬ 
tion of the secondary stars in 2A1822-371, MXB1636-536, MXB1735-444 
and XTE J1814-338 and the determination of their orbital velocities, which lie 
in the range 200-300 km s“^ (see Casares et al. 2003, Casares et al. 2004b). 
In addition, the application of this technique to the BH candidate GX339-4 
during its 2002 outburst provided the first determination of its mass function 
and hence dynamical proof that it is a BH (Hynes et al. 2003). 



Time (ms) 

Figure 5. Millisecond pulses detected 
during the rise of an X-ray burst in 
MXB1636-536 (after Strohmayer et al. 
1998). 



Figure 6. Time evolution of the 580 Hz 
oscillation in MXB 1636-536 during a su¬ 
perburst episode. The solid curve shows 
the best orbital model fit which constrains 

the NS velocity amplitude to 90-175 km 
^-1 


LMXBs are considered to be the progenitors of BMPs because they provide 
a mechanism for spinning up NS to millisecond periods, through the sustained 
accretion of matter with high angular momentum during their long active lives. 
Despite intensive efforts over 2 decades, the detection of millisecond pulsations 
in LMXBs proved elusive but the advent of XTE changed things dramatically 
with the discovery of: (i) persistent pulses in 5 transient EMXBs with spin 
periods in the range 185-435 Hz (ii) nearly coherent oscillations during X-ray 
bursts in 13 EMXBs. Eigure 5 presents an example of a train of oscillations, 
with a frequency of 580 Hz, detected during an X-ray burst in MXB1636-536. 
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In two LMXBs, SAX J1808-3658 (Chakrabarty et al. 2003) and XTE J1814- 
338 (Strohmayer et al. 2003), burst oscillations were detected in addition to 
persistent pulses and they showed identical frequencies. This confirmed that 
burst oscillations are indeed modulated with the spin of the NS. 

Moreover, burst oscillations were detected in an 800 s interval during a su¬ 
perburst in MXB1636-536. The oscillations showed a clear frequency drift 
which was attributed to the orbital Doppler shift of the NS (see fig. 6). A cir¬ 
cular orbit model was fitted to the data and used to constrain the projected NS 
velocity to between 90 and 175 km s“^ (Strohmayer & Markwardt 2002). This 
is a remarkable result which shows that burst oscillations can be used to trace 
NS orbits and, in combination with information provided by the Bowen fluo¬ 
rescent lines, turn persistent LMXBs into double-lined spectroscopic binaries. 

3. Conclusions 

In the past 15 years the field of X-ray binaries has experienced signifi¬ 
cant progress with the discovery of 16 new BHs and 6 millisecond pulsars 
in LMXBs. Reliable mass determinations have been provided for 15 BHs and 
33 NS which are starting to reveal the mass spectrum of compact remnants. 
Two X-ray binaries, V404 Cyg and GRS 1915-1-105, contain BHs too massive 
to be explained by current theoretical modelling. However, better statistics 
are required in order to use the observed distribution to constrain fundamen¬ 
tal parameters of X-ray binary evolution and supernova models, such as the 
mass-loss rate in the W-R phase, the explosion energy dependence on progen¬ 
itor mass, the amount of fall back or details of the Commom Envelope phase. 
The discovery of fluorescence emission from the companion star, together with 
X-ray burst oscillations has opened the door to derive NS masses in persistent 
X-ray binaries. This new technique, which will benefit from new instrumen¬ 
tation on large telescopes (e.g. OSIRIS on GTC), will likely provide further 
evidence for the existence of massive NS in LMXBs. 

Note added. Two new NS spin periods were discovered when finishing 
writing this contribution: burst oscillations at 45 Hz in the persistent LMXB 
EXO 0748-676 (Villarreal & Strohmayer 2004) and persistent pulsations at 
598.9 Hz in the new transient LMXB IGR J00291-I-5934 (Markwardt, Swank 
& Strohmayer 2004), discovered by INTEGRAL in December 2004. 
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